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INTRODUCTION 


This document is the Summary Progress Report submitted by Robert Gold 
Associates in accordance with the CDRL of Contract No. N00173-80-C-0480 and 
describes the work accomplished over the reporting period of September 17, 1980 
through July 1, 1981 on the Study of Multistate Pseudonoise Sequences. 

Although mainly binary encoding sequences have been used to date as 
the encoding mechanism for spread spectrum systems, and their employment has 
resulted in striking improvements in performance and operational capabilities 
over conventional communication systems, it has been recognized that further 
substantial improvements in spread spectrum communication systems could be 
economically achieved through an increased knowledge and application of multi¬ 
state (as distinguished from binary) sequences. 

- The present work has been successful in extending the fundamental con¬ 
cepts and results of the theory of binary sequences to the multistate case. 
Substantial progress and new results have been achieved in the following areas: 

1. Generation of Maximal Linear Multistate Sequences 

2. Implementation of Multistate Sequences with Binary Logic 
Elements 

3. Special Properties of Multistate Sequences 

4. Definition and Computation of the Autocorrelation Function 
of Multistate Sequences 

5. Computation of the Power Spectral Density of Maximal Linear 
Multistate Sequences • 

In Section 1 the multistate encoded signal and its analytical repre¬ 
sentation is reviewed to ensure that the properties of the multistate sequences 
which have the greatest impact on the resultant signals were emphasized during 
the current study. 

In Section 2 the concept of a linear multistate sequence is discussed. 
The multistate alphabet (finite fields over (IF(2 )) and the polynomial corres¬ 
ponding to a feedback shift register are described. A method for implementing 
multistate sequences with binary logic elements is presented in this section. 
This result assures the practical and economical application of multistate 
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sequences. The concept of a maximal multistate sequence is discussed, and a 
method for determining maximal multistate polynomials is given. A table of 
irreducible polynomials of degree 2 for the 8-state case is constructed and 
presented in Section 2.5. 

The autocorrelation function of multistate sequences is defined in 
Section 3, and the autocorrelation function of maximal linear multistate 
sequences is determined. The result is essentially independent of the period 
of the sequence and depends only on the number of states of the sequence. It 
is found that the autocorrelation function of multistate sequences has regu¬ 
larly spaced, relatively'high out-of-phase autocorrelation peaks. In Sec¬ 
tion 3.8 we describe how these out-of-phase correlation sidelobes of a multi¬ 
state maximal linear sequence may be used as an aid to rapid synchronization 
of pseudonoise encoded systems. 

The power spectral density of maximal linear multistate sequences is 
determined in Section 3.8 where a formula for the spectral lines is given in 
terms of the discrete Fourier transform of a sequence of 2 - 1 values, where 
2^ is the number of states. This formula allows the easy computation of the 
spectrum of multistate encoded signals. The power spectral density of a maxi¬ 
mal 8-state sequence of period 63 is computed and presented in this section. 


2 








1. DEFINITION OF MOLT 1 STATE SIGNAL AND ITS CORRELATION FUNCTION 

Although tlie current study is concerned with the properties of the 
digital multistate sequences, which are the Sasic encoding mechanism for 
spread spectrum signals, it is important to define the role of these multi¬ 
state sequences in determining the properties of the resultant signals. This 
will ensure that the direction of the current work can give priority to the 
study of those properties of the multistate sequences which have the greatest 
impact on the resultant signals and their application to tactical communica¬ 
tion systems. Accordingly, a review of the analytical representation of 
pseudonoise signals and their correlation function was undertaken and docu¬ 
mented. The basic formulas are summarized in what follows. 


A pseudonoisc signal which has been encoded by a 
multi state pseudorandom sequence may be represented 
analytically as 

iV* 

s(t) = 2 COS (lot + Y(a(k))) f 0 (t-Ak) 

k = - » 


Re 


i wt 


2 e 

k= - oo 


1 yC 3 ( k ) 3 f Q ( t _ Ak) 


where m = carrier frequence 

a = multistate sequence 

Y = real valued mapping of multistatc sequence 
which determines how the multistate sequence 
will operate on the carrier 

f = squarewave function; f () ft) = 1; 0 < t < A 

f () ( t 1 = 0 ; otherwise 

A = code chip width. 
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Example: Let a be the 4-state sequence with values 

0, 1, 6, B 2 . 

Let y be the mapping such that 
Y : 


B 


0 

1 

B 

2 


7T 

4 


- t r 

X 

3tt 

4 


IT 


The resultant phase modulated signal is illustrated 
in Figure 1. 



Figure i. Four-phase modulated carrier 

The correlation function of wideband PN signals is 
a crucial parameter in measuring their performance in 
spread spectrum systems. The cross-corre 1 ation function 


















of two pseudonoisc signals s,(t) and s,(t) which have been 

n d 

encoded by two multistate sequences a and b, respectively, 
has been shown to be given by 

oo j . P " 1 *_ 

R(s a ,s b )(T) = ^ 4Re|c la,T i^^n(a(k))n(b(k-il))jR(f 0 ,f 0 )CT-Ak) 

where 

n(a)(k) = and n(b) = 

in view of this result, the correlation function of 
two multistate sequences has been defined as follows [1, 
p. 97]. 

Definition : Let a and b be two multistate sequences 

of period p. Let n be any complex valued mapping of the 
sequence alphabet. The correlation function of the multi- 
state sequence with respect to the mapping n is defined 
to be 

t P-1 - 

0 n ( a , b ) (r) = i. y n (a(k) ) n (b (k- t) } • 
n p k = 0 

Hi us , 

CO 

R(s a ,s b ] = y i-Re{o lwT 0 (a,h)(e)}R(f ,f )(r-Al) 

I = - “ * 


1 - 1 Detection of Signal Encoded with Multistate Signal 

The received signal is given by 


s 

a 


(t) 


E cos/w^t + Y(a(b)) 
k=-^> V 


f 0 (t - Ak) 



i w t ,r ' 

c E 


iV(a(k)) 


1 


e 


e 


• f^it - Ak) 










Let the locally generated signal be 


S b ( t-T) = 


i cos(( 


(^c — oj i) t + y (b (k))) f n (t-x 


>) V 


where w ^ is an IF frequency 


y Re 


i (w - w,) t 
e L e 

k=-°° 


iY(b(k)) 


f o (t 


We first multiply the received signal s (t) by 
generated signal s ^(t- - t) 


s (O 
a 


>- 


s (t) s (t 
a b 


s b (t- T ) 


Le t 


iY (a(k)) 

U a (t ) = I) e f n (t - Ak) 

k=- w 


0 


U fe (t) = 12 

k=-°° 


iy(b(k)) 


f Q (t - Ak) 


Then we have 


s a ( t > = 2 R ° 


iw t 

e C U (t) 
a 


s b (t) 


Re 


i(w - at ) t 
e C U b (t) 


- Ak) 

t - Ak) 

the locally 

- t) 


s (t)s,(t-x) 


. iw t i(ui - u>. ) t 

7 Re e C U (t) Re e C 1 


U, (t - t) 











Using the relationship 


Re ( z 1 ) Re ( z 2 ) 




Z 2> 


+ Y Re ( z ! ' V 


where z^ and z 2 are complex numbers and bar indicates complex 
conjugation, we have 


S a {t)s b (t - T > 


8 RC 


i ( 2u) - to i ) t 

c i 


U ( t) U (t- T ) 
a b 


+ -Re 


U (t)U, (t- T ) 
a o 


We now bandpass filter about the IF frequency u) 


to obtain 


y (t) 


j_ 

8 


Re 


i u 
e 


t 

U ( t)U. ( t — T ) 
a b 


and bandpass liLtor y ( t) witii a filter whose lowpass equivalent 
computes 


AP 

/ U a (t)U b (t- T )dt 

0 


Ibis integral has been shown to be equal to 


T: e(a,b)(i) 0 ( f 0 , f 0 ) (t-AI) 

J — — m 

For x = ilk, the above sum becomes 0(a,b)(k) and the output of 
the last bandpass filter is 


7 ( t) 


8 


R e 



t 

9(a,b)(k) 


To envelope detect the waveform z(t) we square to obtain 


, 1 
128 , 

Ro J 

i 2u t 2 

e 9 ( a , b ) ( k ) 

+ Re 

0 ( a , b ) ( k ) 

2~] 

< 

i 1 

J 





and lowpass f i 1 t e r to obtain -1 a - > ^J -L H ) .. 



















































































































































2.2 Multistatc Feedback Shift Registers 

Multistate sequences generated hv feedback shift 
registers may be represented by polynomials with coeffi¬ 
cients in the alphabet field of the sequence. Examples 
are given below for the binary, 4-state, and 8-state cases 

B i n a r y 















8 - S t a t e 



AAAAA 


a 


1 

f(x) 


6 

2 

1 + x + gx" 


= A + /b\ x + Ax 2 + Ax 3 + 

The polynomial corresponding to a multistate shift 
register determines the properties of the multistate 
sequence generated by that shift register. Thus, to find 
the period of any polynomial f over G! ; (2 n ), multiply the 
polynomial by its conjugate polynomials to obtain a poly¬ 
nomial f' over GF(2). The period of f is equal to the 
period of f'. 

bxamp 1 e_ ( 4 - s t a t c) : Gbf2 2 ) f(x) = x“ + ax + 1 

2 7 

The conjugate polynomial is x + a x + 1 


7 2 2 4 12 

(x“ + ax +l)(x + a x + 1) = x +x + x + 


4 .1 

I lie period o I x + x + 
■> ? 
of x" + ax + 1 (and x “ + 


1 

x" + x + 1 is 5 and hence the period 

? 

a"x + lj is 5, e.g., 


T 
X “ 


+ 


The 


ax + 1 


bxamp l e 


con i tig a t e 


2 3 S 0 

+ ax + ax*” + x' + x' + ax + = laalO | laalO | 

( 8 - s t a t c) : Gb (2 3 ) f (x) = x 2 + a 5 x + 1 

2 3 2 6 

po 1 vnom i a 1 s are (x“ + a x + 1) and (x + a x + 1) . 
+ 1 ) ( x “ * a' 1 x + 1 j (X “ + a ( ’ x + 1 ) = x h + X 3 + 1 . 
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The period of x + x' + 1 is 9 and hence the period of 
? 5 

f(x) = x“ + u x + 1 (and its conjugates) is also 9, e.g., 
1 

2 7 5 7 , 
x + a x + 1 

= | 1 a 5 a 2 « a a 2 a 5 1 0 11 a 5 
U-period 9-fcj 


5 f 2 7 3 4 ? s 56 

= 1 + a x + a x“ + ax +ax + a x + a x 


+ x 7 + x 9 + a 5 x 10 


2.3 Polynomials Corrospundito M n I L i s t a t e Shift Registers 

For the case of linear binary sequences, extensive 
tables of irreducible polynomials have been compiled [2] 
which determine the feedback taps for the shift registers 
which generate the sequencer.. 1'he properties of the 
sequences generated by the shift registers may be deter¬ 
mined from the algebraic properties of the corresponding 
polynomial. A method of determining the polynomials for 
multistatc sequences having prescribed properties has 
been developed and is described in this section. The 
method assumes that tables of polynomials for the binary 
case arc available. 

Rcsu1t : To find a polynomial f(x) of degree k 
corresponding to a k-stage shift register which will 
generate a 2 m -state sequence of period p, we proceed 
as follows: 

(1) Assume a is a root of a maximal binary polynomial 
g(x) of degree k*m. 

(2) Choose an integer j such that 

7 m • k j 

^ ~~ rr ->m • k 77 

ged(j ,2 -1) 

k - 1 i *m 

(3) f(x) = TT (x + a J ). 

i = 0 
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Hxamp 1c (4-state): in = 2 


Wo construct a polynomial of degree 2 corresponding 
to a 2-stage shift register which will generate a 2“ = 4 
state sequence of period p= 5. 


(1) Let u he a root of the maximal binary polynomial 

4 4 

of degree m k = 4 given by g(x) = x + x + 1, l. e . , a = a + 1 . 


(2) Choose j = 3 such that 


,in*k ^ -74 _ j 

gcd( j ,2 m ' k - 1) ged (3,2 4 - l) 


(3) f (x ) = (x*a^)(x + a^) = x‘'+(a 3 + a 1 “)x+l 


-7 7 1 |Q 

= x“ + («“ + <* + l)x+l = x" + a x + 1 


7 T $ 

1 (x) = x“ + 6”x + 1, where b = a‘ . 


The shift register corresponding to this polynomial 



The sequences generated by this shift register have 
period 3, c.g., 


I 1 1 6 0 61 1 1 
| period 5 | 

2.4 implementation o f Multistate Sequences 

[n this section, we explain how multistate sequences 
may be implemented with binary sequences generators. We 
shall use an 8-state sequence generator to illustrate the 
general procedure. We consider the maximal polynomial 
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over GF(8) given by f(x) = ax + x + l. The shift register 
corresponding to this polynomial is 



1 1 


2 2 

and the sequence l/(a x +x+ 1) which is generated by 
this shift register with the indicated initial conditions 
is given by 


ax” + x + 1 

The complete sequence and its binary representation is 
shown in Figure 2. 

In order to find the representation for the 8-state 
sequence as three binary sequences, we proceed with the 
following steps. 

Step 1 : Rationalize the polynomial representation 
of the 8-state sequence, e.g., 

1 _ 1 (a^x“ + x+l)(oi^x“ + x + l) 

2 2 * 7 2 4 1 

ax 4 + x + 1 ax + x + 1 (a“x + x + 1 )(a x“ + x + 1) 

6 4 S3, 2 , , 

ax +ax +ax +1 

6 4 3 . 

x + x + x + x + 1 


The denominator of the rationalized fraction provides the 
feedback taps for the required binary shift register. 






r 



Figure 2. Mnxim.il 8-State Sequence and Its Binary Representation 












Stop 2 : Express the coefficients of the numerator 
from GF(8) as binary triples, e.g., 

a 6 = [0 1 1 ] ct = [010] 

a 5 = [1 1 0] 1 = [1 0 0 ] 

a 6 x 4 + '.x 5 x J + ax 2 + 1 = [Pll]x 4 + [110]x J + [0 10 ]x 2 + [ 100 ] 

6 3 3 ” , 6 3 3 , 

x+x+x+x+1 x+x+x+x+1 


Step 3 : Write the fractional representation of the 
sequence given in Step 2 ns the sum of three binary frac 
tions, e.g., 

[011 ]x 4 + [1 10]x ^ + [( ) 10 ]x 2 + [100] 

0 4 3 . 

x + x + x + x + 1 


x~* + 1_ 

x b + X 4 + X 3 + X + 1 


= 1111110 1 


4 3 

X + X + X 


0 4 3 , 

X + X + X + X + 1 


= 00101010 


0 4 3 . 

x + x + x + x+ l 


= 0 0 0 0 1 1 1 0 


. 5 , 4 3 6 , 

11a laaa 1 


The implementation of this 8-state sequence generator 
using a simple binary shift register is illustrated in 








Figure 3. Implementation of an 8-statc generator with 
a single binary feedback shift register 

2.5 Construction of Tables of Irreducible Polynomials of 
Degree 2 over G F ( 2 '■*) _ 

In order to illustrate the computation of useful tables 

of multistate irreducible polynomials, and in order to have at 

hand a convenient list of such polynomials to aid in tiie study 

of multistate sequences, a table of irreducible polynomials 

3 

of degree 2 over the finite field GF(2 ) was constructed. 

These polynomials will generate eight (8) state pseudonoise 
sequences whose periods are factors of 63. 

The multistate polynomials in each group of three (3) of 
Table 1 are factors of the binary polynomial in the product 
column. T hu s, for example 

(x 2 + p. 3 x + B)(x 2 + p/’x + P 2 )(x" + P 3 x + B 2 ) = x 6 + x + 1 

16 










3 

Table 1. Irreducible Polynomials of Degree 2 over GF(2 ) 


Period 


Product 



2 


3 

1 

X 

+ 

6 x + P 


2 


6 

2 

X 

+ 

6 x + P 


2 


5 

4 

X 

+ 

6 x + 6 


2 


5 

3 

X 

+ 

e x + e 


2 


3 

6 

X 

+ 

6 x + 6 

12 

2 

X 

+ 

P 6 x + 6 


2 


2 

5 

X 

+ 

ex + e 

10 

2 

X 

+ 

e 4 x + e 

20 

2 

x 

+ 

Px + e 6 


7 

14 

28 


5 

x + e x + l 


X 4 + ? X + 1 


x 2 + £ ^ x + 1 


63 


21 


63 


x + x + 1 


6 ^ 4 2 , . . 

X + X + X + X+1 


6 , 5 , 2 , , , 

X + X + X + X+1 


X 6 + X 3 + 1 


11 

22 

44 


x 2 + x + g 4 


x + x + 6 

2 2 
x + x + f 


63 


6 5 3, 2 . 

x + x + x + x +1 


1 3 

26 

52 


x 2 + e 6 x + e 6 
x 2 + e 5 x + e 3 
x 2 + e 3 x + e 3 


63 


6 . 5 3 2 . 

X + X +x + x +1 


15 

30 

60 


x + 6 x + B 

2 2 

x + ex + g 

2 2 4 

x^ + e x + e 


21 


6 A 5 4 2 . 

X + x + x + x +1 


21 


X + X + 1 


X + X + 1 
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Table 1. Irreducible Polynomials of Degree 


I 



2 


4 ' 

23 

X 

+ 

6 x + 6 

46 

2 

X 

+ 

4 

6x + 6 


2 


.2 

29 

X 

+ 

6 x + 


2 


„ . 3 

31 

X 

+ 

6x + 6 


2 



62 

X 

+ 

6“x + 6 

61 

2 

X 

+ 

4 

e x + e 


63 


63 
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The roots of the polynomial x + x + 1 are 

, 2 A 8 16 32, 

< a a a a a a J 


and a 

8 



of 

2 


3 




are 

root s 

x 

+ 

£ 

+ 

B 


and 

16 



of 

2 


6 


„ 2 

9 

a 

are 

roots 

X 

+ 

B 

+ 

C 

B = a 

and 

32 



of 

2 


5 


U 

3 „ 2 

a 

are 

roots 

X 

+ 

£ 

+ 

a 

£ + £ + 1 = 0 


2.6 Some Properties of Multistate Sequences 

2.6.1 Multiplication of a Multistate Sequence by a Character 

The following result describes what happens when a multi-state sequence 
is multiplied by a character. The resulting sequence is a phase shift of 
the original sequence. The precise phase shift obtained is described by 
the following result. 


Result: Suppose we have a 2 -state maximal sequence h corresponding to 

the n th degree polynomial f(x) e GF(2 k )[x] with root a V in GF(2 kn ). 

2 k - 1 

(The values taken by this sequence are {0, 1, 6 ... £ ) 

2 k ‘ n - l M i -1 

where £ = a —-- = a .) Then 8 J h (i) = h(i - -jv M) . 

2-1 


Proof: By previous result there exists a linear mapping 

k k ~i 

if of GF(2 -n ) into GF(2 ) such that h(i) = f((a % ) ). 


_ 

1 

- jv 

Mj = 

-(i 

- jv 

a V ) 

-i *1 

a V ) 

, Ml 

(a )! 

-i 

n 

a V ) 

£ J = 


= £ J h (i) 


Example: (see following page) 
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2.6.2 Binary Representation of Maximal Multistatc Sequences 

3 

Result: Suppose we have a 2 -state maximal sequence h corresponding to 

the degree polynomial f(x) c GF(2^)[x] with root a V in GF(2 3n ). 
g (x) 

Let be a rationalized representation of h. Then the binary 

representation of the sequence h = is given by 


h = M- p n , * p i h 

' 7 7 — (v M) ' 7 


(v M) 


where: M = 


2 3n - 1 
2 3 - 1 


Proof: h = 


'Pj(g(x)) P 2 (g(x)) P 3 (g(x))] 


[Pjigfx)) Pj (£ _1 (g(x) )) PjC6(g(x))) 


0} p i( (e_1> (ffS))- 


p . ,f) frS 




* b )-( F ' h ),- m - ( p > h ) 

-(v M) 


Example: Suppose we consider the 2 -state maximal sequence h correspond¬ 

ing to the 2 n3 degree polynomial x" + a^x 4 a^. The sequence gener¬ 
ated by the polynomial and the binary representation of the sequence 
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have been shown in Section 2.4 to be the sequences of Figure 2. We note 
from Table 1 that a* 3 is a root of the polynomial x 2 + a^x *- . Thus 

v = 13; v =34 and v 3 M = 54. Thus the binary sequence representation 
of the 8-state sequence is 



This may be verified by observing the binary sequence representation of 
Figure 2. 

2.6.3 Sequences Generated by Conjugate Polynomials 

R 

The k-1 conjugate polynomial to a polynomial f(x) e GF(2 ) are obtained 

2 

by replacing the coefficient f(j) of x^ by (f (j )) *" , (f(j)) 2 , . . . 

2 k_1 

f(j) *■ , respectively. Thus, for example, the conjugate polynomials to 

x 2 + B^x + B^ c GF(2 3 ) [x] are found to be 


l 2 x 2 + (B 6 ) x + ( f V° - x 2 + B 5 x + B 5 


l 4 X 2 + (fS x + (f 6 ) = X 2 + e 3 + B 3 


The 2 state sequences generated by these conjugate polynomials can be 

obtained from the sequences h generated by the polynomial f(x) by replacing 

2 k-1 

each term h(i) of the sequence by (h(i))~, (h(i))~ , . . . (h(i)) 2 , 

respectively. Thus in the above example if h(0), h(l), ... is the 

sequence generated by f(x) = x“ + u^x + , then (h(0))", (h (1)) . . . and 

(h(0)) 4 , (lid)) 4 , . . . will be the sequences generated bv x~ + a 3 x + a 3 and 
o 3 3 

x + a x + n , respectively. 
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Since the j 1 * 1 conjugate polynomial to the polynomial f(>:) f C.F(2 ) 

, j t 

having root a X has root v , the binary representation of the sequence cor¬ 
responding to this conjugate polynomial is 



^66 13 

Example: The polynomial x + £ x + £ having root a will generate the 

5 A 3 

sequence h = 1 1 £ 1 £ £ £ .... This sequence can be repre¬ 

sented by phase shifts of a binary sequence a corresponding to the 


polynomial x* 3 + x - " + x 3 + x 2 + 1 = (x 2 + £ x+£ ) (x + £ 3 x+B ) (x +B x+B 3 ). 


In fact in Figure 2 we have 


h 

= i 

1 

„5 

D 

i e A 

e 3 

a 

= i 

1 

1 

i i 

i 

a 9 

= 0 

0 

1 

0 1 

0 

a „ 

= 0 

0 

0 

0 1 

i 

-y 

The conjugate 

polynomial 

x 2 + £ 5 x 

f £ 6 

sequence 

1 

h = 

1 1 

B 6 1 F 

£ 6 £ 5 . 

• • 9 

sequence 

h by 

replacing 

each term h(i 


0 ... 

1 ... 

1 ... 

2 6 

having root a will generate the 
which is obtained from the 
by (Vi(i)j". The sequence h 
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can be represented by phase shifts of the sequence a in accordance with 
the above result. Thus we have 



? 

h 

1 

1 f 3 

1 

e e 6 

E 5 

. . . 

a 18 

t 

a = 

1 

1 1 

1 

0 0 

1 

. . . 

a 32x9 

a 3fc = 

0 

0 0 

0 

l i 

1 

. . . 

a -32x9 

a -36 = 

0 

0 1 

0 

0 1 

0 

. . . 

The conjugate 

polynomial 

2 

X 

3 3 

+ ex + e 

having root 

52 

a ’ 

will generate 

the sequence h 

t 1 

= 116 

6 1 

2 5 3 

e e b . 

• • 

. which is 

obtained from the 

sequence h by 

replacing each 

term h(i) 

by 

(h(i))\ 

The 

1 F 

sequence h 

can be represented by phase 

shifts of 

the 

sequence 

a in 

accordance 

with the above 

result. 

Thus 

we have 






t t 

h 

1 

i e 6 

1 

b 2 e 5 

e 3 

. . . 

a 27 = 

1 f 

a = 

1 

1 0 

1 

0 1 

i 

. . . 

a 16x9 = 

t i 

a l 8 

0 

0 1 

0 

0 1 

0 

. . . 

a -16x9 = 

« f 

a -l 8 = 

0 

0 1 

0 

1 0 

i 

. . . 


v 


♦ 
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3. 


CORRELATION FUNCTION OF MULTI STATE SEQUENCES 


V 

i 


I 


3.1 Definition of 9 

The correlation function with respect to a complex valued mapping of the 
sequence alphabet of two (2) multi-state sequences has been defined (in Report 
No. 1) as follows 


6 (a ,b) (t) 
n 


t’-i 

= ^ ^ ^ n(a(k)) r(b(k-r)) 


where a,b are multi-state sequences of period P 

n is complex valued mapping of the sequence alphabet 

In what follows, we compute the auto correlation function with respect to 
a specific complex valued function n of maximal multi-state PN sequences. 

3.2 The Mapping of n 

The specific complex valued mapping of CF(2 ) which we shall use has 
values on the complex unit circle. The mapping n ultimately determines the 
correspondence between the characters of the multi-state alphabet and the phase 
of the carrier which will be modulated by the multi-state sequence. 

Definition: n is a complex valued mapping of 

CF(2 k ) = jo, 1, B, E" . . . 6 2 “ > given by 


n (0) 1 


n(6 s ) 


We note that 


i (s + 1 )tt 
. k-1 ' 


0 < s < 2 -2 


n(B S ) 


i (s + 2)m 
~k-1 


2 k - 1 1 s i 2 k+1 - 3 











■? fc ' 









6(a,a) (t) 


-1 


otherwise 


2 n ‘ k -l 


6 (.a,a) (t) = n(a(t)^ n(a(t-t)j 

t=0 


definition of £ 


6 (a,a)(jM) 


2 n ' k -1 


^ ^ n (a(t)) n(a(t-jM)J 


t=0 


2 n ' k -1 


^ ^ n (a (t)) 


n a (t - (jv)(v 


■'><)) 


t=0 

2°' k -1 


^ ^ n(a(t)) a(t)^ 


t=0 


M = 


2 k n -I 
2 k -l 


previous result /' (3.0) 


For each non-zero term a(t) e GF(2 ) of the multi-state sequence a 
there is an integer f ^a (t)) where 0 2 E (a (t jj - 2 k -2. Such that 


(a(tl) 


a(t) = 6 V > . Thus we have 


2 n " k -l 


y^n(s Ua(t)) )r,(e jv + t(a(t)) ) + 2 k 


(n-1) 


- 1 


t=0 

a(t)/0 


2 n-k -l 

£ 

a(t)*0 


i(t(a(t)) + l) it 

7k-1 


i (f (a (t)) + v j + 7i 


,k-l 


jv+ JL(a (t) ) < 2 -2 
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2 nk -l 

e 

a(t)^0 


i 

J(a(t)) + l)* 

i (i (a (t)) + vj + 2 y 

71 

e 

2 k-1 

2 k-l 

e 



+ 2 k(n-l) _ j 


jv+ £.(a(t)) > 2 -2 


2 n ’ k - 1 

■ E 

a(t)/0 
i(a(t)) 1 2 k —2—j 


i v j it 

2 " -1 

V' 

i (vj + 1 )ir 

2 k_1 

+ X 

2 k_1 

a 

o 

II 

<3 

e 

u 


+ 2 k(n - 1} - 1 


£(a (t)) > 2 -2-vj 


k k 

Since there are (2 -l)-vj values for which J(a(t)) i 2 -2-vj, and 

there are vj values for which £(a(t))> 2 -2-vj, and since each value 

.k(n-l) . 

occurs in the sequence 2 times, we have 


0(a,a)(jM) 


( 2k -‘-vj) 


k(n _ I} _i^p k(n . 1} . JL(vj±i)l k(n _i) 

2 i / ., 2 


+ (vj) 2 e 


+ 2 


3.4 Autocorrelation Function for 8-State Maximal Sequences 

For 8-state maximal sequences, the above correlation formula becomes 


2 3n -l, 


2 3 -l 

©|.i J= (7- vj) 2 e 


3(n-J) 3 (n — 1) . iMtllZ. 3 (n-l) 

A M 


+ (vj) 2 e 


+ 2 


1 


For a two-stage maximal register which generates an eight-state sequence 
of period 63 (n=2), the correlation formula becomes 


i v j t i (vj + 1 )tt 

4 4 

6(9j) = (7-vj)8e + vj 8 e +7 


Example: We use the above formula to compute the auto correlation function 

of the maximal 8-state pseudo noise sequence corresponding to the 
polynomial x 3 + a^x + a 3 . We note from Table 2-1 that a 33 is a root 
of this polynomial. We thus have 
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i 2 i ti 



(21+ l)w 

e(9j) = 

(7 - 2j) (8) e 4 + 

(2j)(8) e 

4 + 7 




1271 

1 3t 



e (9) = 

(5)(8) 

e 

4 + (16) e 

4 + 

7 = 

(7 - 8/2) + i(-40 - 8/2) 




i4r 

i5*n 



6(18) = 

(3)(8) 

e 

4 + (32) e 

4 + 

7 == 

(-17 - 16/2) + i 1 bfi 




i6r 

i 7 TT 



e(27) = 

(1)(8) 

e 

4 + (6)(8) e 

4 + 

7 = 

(7 + 24> / 2) + i (8 + 24.2) 




ir 

i 2rr 



6(36) = 

(6)(8) 

e 

4 + (1)(8) e 

4 + 

7 = 

(7 + 24/2) + i (-8 - 24/2) 




1 3tt 

i4r 



6(45) = 

(4)(8) 

e 

4 + (3)(8) e 

4 + 

7 = 

(-17 - lb.'T) + i(-16. / 2) 




i5r 

6 r 



6(54) = 

(2)(8) 

e 

4 + (5)(8) e 

4 + 

7 = 

(7 - 8.2) + i (40 + 8.1) 


Example: We use the above formula to compute the auto correlation function 

of the maximal 8-state pseudo noise sequence corresponding to the 

2 4 2 23 

polynomial x + 6 x + £ . We note from Table 1 that a is a root 

of this polynomial. We thus have 

_ i2 i~ _ i (2.i + 1 )t 


e (9 j) 

e (9) 

6(18) 

6(27) 


(7 - 2j)(8) e + (2j)(8) e 

i 2 it i 3 tt 

(5) (8) e 4 + 16 e 4 + 7 = 

i4r i5r 

(3)(8) e 4 + 32 e 4 + 7 = 

ifer i 7 tt 

(1)(8) e 4 + (6)(8) e 4 + 7 = 


+ 7 

(7 - 8 ^) -t i (-40 - 8/2) = (-4 • 31 - i (51 • 31 j) 
(-17 -16» / 2‘) + 1(16/2 = (-39-63+i 22 • 63) 


(7 + 24*'2) + i(8 + 24/2) *= (40 • 94 + i 41-94) 








| 0(a,a)| 



Figure U. Absolute Value of Autocorrelation Function of Maximal 
Linear Sequence of Period 63 over GF(2 3 ) 


Sequence Period 



Autocorrelation Function of Maximal Linear Sequence over GF(2 ) (Eight State) 















XTT 12^ 

6(36) = (6) (8) e 4 + (1)(8) e 4 + 7 

i 3tt i4r 

6(45) = (A)(8) e 4 + (3)(8) e 4 + 7 

_ i5tt i6r 

6(54) = (2)(8) e 4 + (5)(8) e 4 + 7 

A plot of the absolute value |6(a,a)| is presented in Figure 4, and of 
9(a,a) in Figure 5. 

3.5 Computation of the Autocorrelation Function of 16-State Maximal Linear 
Sequences _ 

In this section we give examples of the computation of the autocorrela¬ 
tion function of a 16-state maximal linear sequence. The unnormalized auto- 

k k • n 

correlation function of a 2 -state maximal linear sequence of period 2 - 1 

has been shown to be given by the formula: 

0(a,a)(jM) = 2 k ^ n '^ Q k (rj) - 1 

9(a,a)( t ) = - 1 r £ jM 


= (7+ 24/2) + i(-8 - 24/2) = (40 • 94 - i 41 • 94) 

= (-17 - 16/2) + i(-16/2) = (-39 • 63 - i 22. • 63) 


= (7 - 8/2) + i (40 + 8/2) = -4 • 31 + i 51 • 31 


where M 


2 k ' n - 1 
2 k - 1 


Sequence Period 


(No. of Sequence States) - 1 


No. of 

Correlation Peaks 


and Q.(j; 


( 2 k - 1 - j) 


i J 

ok - 1 


Ki + U 17 

ok - 1 


+ j e 


+ 1 


Q 4 (j) = 


- 1 J * 

, 8 

(15 - j - l) e + j e 


jliijjz 

8 


+ 1 
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J 


Q 4 (j) 


1 

14.6 

- 

i 

6.1 

2 

11.0 

- 

i 

11.0 

3 

5.6 

- 

i 

14.1 

4 

- 0.5 

- 

i 

14.7 

5 

- 6.4 

- 

i 

12.8 

6 

-10.9 

- 

i 

8.7 

7 

-13.4 

- 

i 

3.1 

8 

-13.4 

+ 

i 

3.1 

9 

-10.9 

+ 

i 

8.7 

10 

- 6.4 

+ 

i 

12.8 

11 

- 0.5 

+ 

i 

14.7 

12 

5.6 

+ 

i 

14.1 

13 

10.9 

+ 

i 

11.0 

14 

14.6 

+ 

i 

6.1 


The above values of determine the auto correlation function of 16-state 
maximal linear sequences of any period. Thus, for example, for n = 2, 
period 255 we have: 


j 

9 (17 J) 


16 q 4 (j) - l 

|6(17j)| 

19(17j)1/255 

1 

232 .6 

- 

i 

97.6 

252.2 

.99 

2 

175.0 

- 

i 

176.0 

248.2 

.97 

3 

88.6 

- 

i 

225 .6 

242.4 

.95 

4 

- 9.0 

- 

i 

235.2 

235.4 

.92 

5 

-103.4 

- 

i 

204.8 

229 4 

.90 

6 

-175.4 

- 

i 

139.2 

224.0 

.88 

7 

-215.4 

- 

i 

49.6 

221.0 

.87 

8 

-215.4 

+ 

i 

49.6 

221.0 

.87 

9 

-175.4 

+ 

i 

139.2 

224.0 

.88 

10 

-103.4 

+ 

i 

204.8 

229.4 

.90 

11 

- 9.0 

+ 

i 

235.2 

235.4 

.92 

12 

88.6 

+ 

i 

225.6 

242.4 

.95 

13 

175.0 

+ 

i 

176.0 

248.2 

.97 

14 

232.6 

+ 

i 

97.6 

252.2 

.99 


The absolute values of 9 are plotted in Figure 6. 
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4 


The expression for the normalized auto correlation function is: 


9(a,a)( jM) = 


2 KM - 

kn 

2 - 1 


This expression may be approximated by: 


9(a,a)(jM) - 


Q K (rj) 


We note that this expression is independent of n. 


3.6 Asymtotic Expression for 9^ 

The absolute value of the normalized auto correlation peaks of multi¬ 
state maximal linear sequences approach 1 as the number of states increases. 
In fact we have: 


Proof: 


9(jM) 

9(0) 


9(jM) 

9(0) 


-> 1 for large k 


q k (j)' 


1 3 T _ M, j + 1 lz 

k \ ? k-l * ? k-l 

2 - 1 - J e + j e + 1 
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1 j TT 1' j TT 

i - 1 


/ k \ J<-1 jc-l 

1 

7T- 

1 


V 2 - V e 2 - j e 2 

1 - e 2 

+ 1 


3.7 Use of Correlation Sidelobes for Signal Acquisition 

In this section we describe how the out-of-phase correlation sidelobes 
of a multi-state maximal linear sequence may be used as a synchronization aid. 
We show that the time required for the synchronization of multi-state sequences 
is significantly less than that required for the synchronization of binary 
maximal linear sequences. 

The conventional synchronization techniques for binary sequences compares 
the received code with successive phase shifts of a locally generated replica 
of the received code until phase synchronism of the two (2) codes is achieved. 
If the codes have period N and the integration is carried out over W code chips 
for each phase comparison, then this synchronization process will require inte¬ 
gration over as many as WN chips. 

The suggested technique for the synchronization of multi-state maximal 
linear codes makes use of their relatively high and regularly positioned auto 
correlation sidelobes. The synchronization procedure uses the conventional 
search technique to locate one of the set of equally spaced correlation side¬ 
lobes. Since the auto correlation sidelobes are spaced M chips apart, at 
most M relative phase shifts of the two (2) codes must be compared in order 


* 


36 












to find one of the auto correlation sidelobes. The location of one of the 
auto correlation peaks determines the phase position of the other regularly 
spaced correlation peaks. Each of these phase positions is then tested to 
determine phase synchronism between the received and locally generated codes. 

Since the correlation sidelobes which are being searched out in this pro¬ 
cedure have a normalized value 8 of less than one, the integration process 
must be taken over a greater number of code chips in order to compensate for 
the fact that the two (2) code phases being compared are not perfectly corre¬ 
lated. Thus if the integration is performed for time T in the case of perfect 

correlation (0 = 1), then the integration time must be T, = in order to 

8 

obtain an equivalent amount of energy out of the correlator. The following 
argument established this fact. 

Suppose that in the case of perfect correlation the voltage out of the 

2 

correlator is V or equivalently the normalized power is V . If we integrate 

p 

over a time T, we obtain an energy level of E = V T. If the normalized corre¬ 
lation function 0 is less than one, then the voltage at the output of the 
correlator is 0V or equivalently the normalized power is (0V) . In order to 

obtain the same signal energy at the output of the correlator as in the pre- 

2 2 

vious case, we must integrate over a time T, such that (0V) Tj = V T or 



We now compare the acquisition times required for binary and multi-state 
maximal linear sequences. We assume that the binary sequence is acquired by 
integrating over W chips in each relative phase position so that the acquisi¬ 
tion process will require integration over as many as WN chips, where W is 
the code period. To acquire a multi-state maximal linear sequence, we first 
find one of the correlation peaks. The maximum number of chips over which 
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we must integrate is chips. The maximum number ot chips over which we 

/ N \ W 

must integrate to find the mean correlation peak is (-w- - 

V M / (1-9) 


1 • 


Thus the 


total maximum number of chips over which we must integrate is 
WM_ ( h \ W 

0 2 \ M / (l-Q) 2 

of chips WN over which we must integrate to acquire a binary maximal linear 


The ratio of this number of chips to the maximum number 


sequence is 


/M \ J_ + J_ 1 
W e 2 " (7^7 


I » fp k * n « 

For maximal 2 k -state sequences N = 2 k n - 1 and M = z — - 


2 - 1 


formula thus becomes 


The above 


a 




2 k n - ijl (i-e) 2 


For 8-state maximal linear sequences we have 


a 



Using the average value of 9 over the six (6) auto correlation sidelobes, we 
have 0 = -83 


a 


•21 + 


243 


,3 -n 


- 1 


This function is plotted in Figure 7. 
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3 

Figure 7 


" I 1 i • 

4 5 6 7 8 

Comparison of Binary and Multistate Synchronization Times 
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3.8 Power Soectral Density of a Maxima] Linear Sequence over GF(2*) 


In this section we determine the power spectral density of a maximal 

k 

linear sequence over GF(2 ) in terms of the complex valued Q-function which 
has been defined as 2 7r i 

Lk , A " j . .. 2 k 


{ k A -j . 

Qj < (j) = (2 - 1 - j+j U) +1 where w = e 

The power spectral density S of a maximal linear sequence over GF(2 k ) will be 


S(o) = 1 

F(Q k )(j) 

S(j) = (P + 1) 

2 k 

0 k 0 2 u i j s 

t v*>« 2l -' 

s=0 

where F(Q k )(j) = 


The form of the autocorrelation function and power spectral density of a maxi 
mal linear sequence is illustrated in Figure 8. 

IS r~ / r\ \ t ' \ 


F(Q k )(j) 


F(Q.) has period 2 K - 1 


2 ( 2 - 1 ) 


3(2-1) 


Figure 8. Power Spectral Density of Maximal Linear Sequence over GF(2 K ) 
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The power spectral density S of a sequence a of period P is found via 
the Weiner Khinchin Theorem by taking the discrete Fourier transform F of 


its autocorrelation function 9 where the Fourier transform is given by 


■ j. 

S(a)(j) = F(9(a))(j) = ]P e(a)(k) e 


2tt i j k 


In this section we show that if a is a maximal linear sequence of period 
P = 2 n ' k - 1, then 


S(a)(o) = 1 


S(a)^j(2 k - 1)^ = 2 kn = P +1 

k(n-l) . 

S(a)(j) = 2 F(Q k ){j) for j is not a multiple of 2-1 

We first make the following computation for Q k 


Result: 


Let Q k (j) = (2 - 1 - j) u> + j w 


Then Q k (j) 
j=0 


•(j+1) 


u = e 


Proof: 


T, * J 5 H 1 - 

j=0 


k k k 

2-1 2 2-1 


X - X - X + 1 
(1 - X) 2 


Differentiating and multiplying by X we have the identity 


k 2 k k 9 k i 


V* j X J _ (2 -2 m - (2 -1)X^ + X 

U U - x) ? 
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Using these identities with X = w we obtain 


2 k -2 

(i) (2 k -1) y 

j=0 


2(2 k -l) - - (2 k -l)ui 

(1 - u )) 2 


2 k -2 

(2) - y j “ j 

j=0 


(2 k -2) -t- (2 k -l)cii~ 1 - u> 

(1 - o >) 2 


(3) 


£ 


j=0 



(2 k -2)u) - (2 k -l) + J 

(1 - u ,) 2 


(1) + (2) + (3) = 1, and hence the result follows. 

We now compute S(a) 

P-1 _ 2rijs 

Proof: S(a)(j) = F(0)(j) = ^ 9(s) e P P = 2 nk - 1 

s=0 


Let 9' be the autocorrelation function of a binary maximal linear sequence 
of period 2 n k -l. It is well known that 


9' (0) 
9*(s) 


2 n " k - 1 F(9 1 )(0) 

and 

-1 for s t 0 F{9')(j) 


1 

P + 1 = 2 kn for j i 
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Since 0(s) = -1 for s, not a multiple of M = 


2 n-k _ j 


2 k - 1 


L , we have 


P -1 2tr f j S 

p 

s=o 


p-i 

2tt i j s 

2-2 

23 e ’^ s ^ e 

p + 

23 D + 0(sM)]e 

s=0 


S=1 

P-1 

2iri js 

2 k -2 

23 8 ’^ s ^ e 

p 

+ 23 [i + e(sM)] 


s=0 s=l 


2tt i j Ms 
P 


2th s 
2 k - 1 


Since 0(sM) = 2 k(n_1 ^ Q k (j) - 1, we have 



•s. 


P-1 

23 0*(s)e 
s=0 


2»ijs 2 k -2 

P + 2 k(n - 1} ^ Q k (s)e 

s=l 


2tt j j s 
2 k -l 


Since Q k (0) = 2 k , 


we have 


P-1 

23 O'OO e 

s=0 


2tt i js 
P 


+ 2 


k(n-l) 


2 k -2 


Z q k (s) 

s=0 


_ 2nijs 
2 k -l 


= F(0')(j) + 2 k(rM) F(Q fc )(j) - 2 k ' n 


i 

- -< 

< ; 
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Since F(9‘)(0) = 1 and F<Q k )(Q) = 2 k , we have 

S(a)(0) = 1 

F(6')(j) = P + 1 = 2 k ' n for j f 0. Hence 

S(a)(j) = 2 k(n_1) F(Q k )(j) for j t 0 

2 k -2 

• Y v s > ■ 2k 

s=0 

Thus S(a)^j(2 k -1)^ = 2 kn = P + 1. 

The power spectral density of a maximal linear sequence of period 63 
is shown in Figure 9. 



44 








10 Log [Power Spectral] of Maximal Sequence over GF(8) Period 63 































APPENDIX A 


Theorem: 


Proof: 


f(O)h(i) 


Let 4> be any linear mapping of the field GF(2°) into the field 

k k 

GF(2 ). Let f be any polynomial with coefficients in GF(2 ), i.e., 

f c GF(2^)[x]. Let a be a root of f in GF(2 n ). The sequence 

h(i) = if (a 1 ) is a member of V(f), i.e., the sequence h can be 

generated by the shif* register corresponding to the polynomial f. 


We assert that the sequence 

h(i) = f(a _i ) 

can be generated by the shift register corresponding to the polyno¬ 
mial f (x) where f(a) = 0. The sequences h generated by the poly¬ 
nomial f(x) are those which satisfy the recursion relation: 

= f (l)h(i-l) + f (2)h(i'-2) + . . . f (n)h(i-n) for all i > n 

Thus, to show that the sequence h can be generated by the polynomial 
f, we must show that is satisfies the recursion relation. 


f(a) = 0 


hence f(0) + f(l)a + f(2)a. 2 + . . . f(n)c. n = 0 


-1 


multiplying by a we have 


f(0)a _i + f(l)a" (l_1) + f(2)cT (1 " 2) + . . . f(n)a" (l “ n) = 0 


♦ (0) 


^f(0)a -i + f(l)a~ (1_1) + . . . f(n)a~ (i ~ n) j = 

(a -1 ) + f(l)f + . . . f(n)^a' (1 " n) ^ - 0 


f(0H 

f(0)h(i) + f(l)h(i-l) + . . . f(n)h(i-n) 


A-l 











Theorem 


k k 

(Converse) Given any 2 -state sequence h with values in CF(2 ) 

which has been generated by a shift register corresponding to the 

irreducible polynomial f(x) of degree n and having coefficients in 

R n k 

GF(2 ), then there is a linear mapping $ of GF(2 ) into GF(2 ) such 

that h(i) * $(a 1 ) where a is a root of f(x). 


A-2 





APPENDIX B 


In this appendix we review the properties of the projection mapping [l]. 

Result: For all Y e GF(2 2 ) we have P 2 (y) - Pj^B *Y^ 

P 3 Cr) = P^BY) 

Proof: Y = a(0) + a(l)B + a(2)B 2 

BY = a(0)B + a(l)B + a(2)(B 2 + 1) 

= a (2) + a (0)B + (a(1) + a(2)j B 2 
6 _1 Y = a(0)B 6 + a(1) + a(2)B 

= a(0)(B + B 2 ) + a(1) + a(2)B 

= a(1) + a(0) + a(2) B + a(0)B 2 

Pj(B _1 Y) = a(1) = P 2 (Y) 

Result: f(x) e GF(2 3 )[x] we have P^ (f (x)) = P j(p *f(x)) 

P 3 (f(x)) = P j (b f (x)) 

Proof: Reference Gold Notes (Page 80) 

P^fU)) 

Pj(b I f (i)x 1 ) 

Pj(l(6f (i)V) 


EP^efU)) x 1 

definition of P 

lP 3 (f(i)) x 1 

previous result 

P 3 (lf(i)x^ 



P 3 (f(x)) 

[l] Study of Multi-State PN Sequences and Their Application to Communica¬ 
tion Systems. Contract N00014-75-C-1040 - 23 July 1976. 


B-l 






Result: 


Let a = 


Proof : 


Result 


Proof : 


g(x) 
f (x) 


be a rationalized representation of a multi-state 


sequence where g(x) e GF(2 )[x] and f(x) e GF(2)[x]. The represen- 

p ( x ) 

tation of a = ■ is the interleaving of k binary sequences is 


given by: 


a = 


P 1 (g(x)) P 2 (g(x)) 


f (x) 


f (x) 


X s(x) ) 

f (x) 


Reference Gold Report (Page 83) 


a 

g 

f 

a • f 
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